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ABSTRACT 


We  consider  the  problem  of  determining  optimal  stock  levels  for  a  multi-item 
inventory  system  in  which  demands  for  items  occur  in  randomly  selected  groups, 
causing  in  interdependent  demand  processes.  This  structure  introduces  dependen¬ 
cies  that  significantly  affect  the  selection  of  item  stock  levels  in  systems  whose 
perfonnance  is  measured  in  terms  of  number  of  jobs  completed  or  time  until  first 
stockout.  Performance  bounds  and  stock  level  selection  techniques  for  several  such 
systems  are  developed  and  analyzed.  ■<= - - 


1  Introduction 


We  consider  a  multi-item  inventory  system  of  the  following  sort.  Demands  arrive 
randomly  in  the  form  of  requests  for  subsets  of  the  multiple  items.  A  demand  is 
met  only  when  its  entire  subset  of  items  is  available.  The  key  difference  between 
this  system  and  many  typical  inventory  models  is  that  it  is  the  number  of  demands 
met  (or  not  met)  that  is  important,  rather  than  the  number  of  items  supplied  (or 
short). 

A  classical  example  of  such  an  inventory  system  arises  in  the  repair  of  machines 
in  the  field,  e.g.,  the  repair  of  computer  terminals,  photocopiers  or  other  equipment 
at  customer  sites.  Service  personnel  are  equipped  with  repair  kits  containing  a 
certain  collection  of  parts  and  tools  that  are  used  for  field  repairs.  Each  repair  job 
requires  some  particular  set  of  parts  and  tools  for  completion.  The  kit  is  typically 
adequate  only  for  some  subset  of  the  possible  repair  jobs.  The  perceived  perfor¬ 
mance  of  the  kit  will  be  the  expected  fraction  of  the  repur  jobs  for  which  it  is 
adequate  or  the  expected  number  of  jobs  which  can  be  completed  before  stockout. 
This  type  of  inventory  system  is  often  called  a  "job  completion”  or  “job  fill”  in¬ 
ventory  system.  Other  job  completion  inventory  systems  vise  in  stocking  spare 
parts  to  assure  the  completion  of  a  specific  mission.  Examples  are  voyages  at  sea, 
space  fiights  or  military  missions.  The  performance  criterion  of  the  inventory  is 
the  probability  that  the  mission  will  be  completed  before  a  required  part  is  not 
avulable  or  the  expected  length  of  time  to  stockout.  Another  application  arises  in 
filling  orders  for  merchandise.  If  orders  occur  in  the  form  of  sets  of  multiple  items, 
the  fraction  of  orders  that  can  be  filled  completely  is  an  important  measure  of  the 
system’s  performance.  An  analogous  problem  arises  in  *build  to  order”  production 
systems  in  which  final  products  have  common  components. 

Job  completion  based  inventory  ssrstems  have  been  considered  in  the  manage¬ 
ment  science  literature  in  a  variety  of  contexts.  Sherbrooke  (1971)  considered  an 
inventory  system  in  which  a  fuled  line  replaceable  unit  (LRU)  could  be  replaced, 
i.e.,  the  repair  job  completed,  only  if  all  failed  subassemblies  were  avulable  as 
spares.  Silver  (1972)  developed  a  dynamic  programming  solution  for  this  prob¬ 
lem.  These  formulations  assumed  that  subassemblies  f^ed  independently  of  one 
another,  according  to  Poisson  procenes.  Smith,  Chambers  and  Shlifer  (1980)  ob- 
tidned  a  closed  form  optimal  stocking  condition  for  this  type  of  problem  in  the 
one  period  case,  if  the  item  stock  levels  are  restricted  to  be  0,  1.  Graves  (1982) 
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reformulated  this  as  a  0,1  knapsack  problem,  which  allows  for  the  inclusion  of  linear 
constraints,  such  as  inventory  space  requirements.  Hausman  (1982)  developed  an 
exchange  curve  analysis  for  this  inventory  stocking  problem  and  Schaefer  (1983) 
integrated  the  independent  demand  job  completion  formulation  into  an  overall  in¬ 
ventory  management  system. 

A  key  assumption  in  the  models  discussed  above  is  that  individual  item  de¬ 
mands  are  independent  of  one  another,  so  that  the  probability  that  any  particular 
combination  of  items  is  required  on  a  job  is  determined  by  the  marginal  demand 
rates  of  the  individual  items.  This  property  does  not  hold  in  many  of  the  examples 
cited  earlier.  In  machine  failures,  for  example,  certain  part  failures  may  contribute 
to  the  failures  of  other  parts,  while  in  order  filling  systems,  demands  may  be  cor¬ 
related  (either  positively  or  negatively)  for  a  variety  of  reasons.  Mamer  and  Smith 
(1982)  developed  an  optimization  algorithm,  using  the  maximum  flow/minimum 
cut  theory,  to  solve  the  0,1  stock  level,  multi-item  case,  when  individual  item  de¬ 
mands  are  arbitrarily  correlated  and  extended  this  approach  to  include  the  use 
of  spare  replacement  machines  [Mamer  and  Smith  (1985).]  March  and  Scudder 
(1984)  developed  an  algorithm  for  determining  optimal  0,1  stocking  policies,  sub¬ 
ject  to  budget  constraints.  Recently,  multi-item  inventories  are  also  being  studied 
in  the  context  of  production  processes  with  parts  commonality  among  products. 
[Baker  (1985),  Baker,  Magazine  and  Nuttle  (1986).]  As  the  examples  analyzed  in 
these  papers  indicate,  determination  of  optimal  stocking  policies  in  such  systems  is 
difficult,  even  for  products  with  only  one  common  part.  Given  the  large  number  of 
part  types  used  in  many  production  applications,  the  ability  to  analyze  larger  prob¬ 
lems  appears  important  in  this  application  area  as  well.  In  a  recent  paper  (1988) 
the  authors  develop  a  lower  bound  heurutic  for  a  model  which  is  q;uite  similar  to 
the  one  considered  here,  but  in  which  the  job  arrival  process  is  a  Poisson  process. 
We  obtain  bounds  on  system  performance  exploiting  the  association  property 
of  part  demand  processes  which  flows  from  the  Poisson  arrivals  assun^tion. 

This  paper  considers  a  formulation  that  is  more  general  than  those  considered 
previously  for  large  problems.  Item  demands  are  dependent  in  an  arbitrary  man¬ 
ner,  as  in  the  more  recent  models  disciused  above.  Stock  levels  are  not  reatricted 
to  0,1  and  multiple  items  of  a  given  type  are  permitted  on  a  particular  job.  Several 
alternative  performance  criteria  are  considered,  including  the  probability  of  com¬ 
pleting  a  fixed  number  of  jobs,  the  expected  number  of  jobs  completed  before  first 
stockout,  the  expected  time  until  stockout  and  the  probability  of  serving  all  jobs 
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arriving  within  a  fixed  period  of  time.  The  resupply  process  for  this  analysis,  as  in 
previous  models,  is  fairly  simple.  It  is  assumed  that  once  a  request  is  unfilled  (a 
job  cannot  be  completed),  all  items  are  restocked  to  their  originally  specified  lev¬ 
els.  This  neglects  the  effect  of  lead  times  and  the  possibility  of  stockouts  at  higher 
inventory  echelons.  This  assumption  appears  to  be  quite  reasonable  for  the  field 
repair  systems  application  discussed  above,  since  repair  personnel  typically  restock 
completely  whenever  they  are  required  to  return  to  the  parts  supply  center.  In  one 
period  inventory  systems  that  support  the  completion  of  a  single  mission,  there  is 
no  restocking  process,  so  this  model  captures  the  relevant  features  of  these  sys¬ 
tems.  For  the  performance  criteria  involving  a  fixed  time  interval,  this  assumption 
is  equivalent  to  a  periodic  resupply  in  which  order  levels  can  be  readjusted  up  to 
the  time  of  delivery. 

This  paper  derives  and  compares  a  variety  of  upper  and  lower  bounds  for  system 
performance.  Since  the  tightness  of  the  bounds  varies  by  application,  a  set  of  upper 
and  lower  bounds  is  advantageous  in  that  the  maximum  lower  bound  and  minimum 
upper  bound  can  then  be  chosen  in  any  specific  example. 

Because  of  the  complexity  of  the  models  considered  in  this  paper,  closed  form 
solutions  for  optimal  inventory  policies  cannot  be  obtained.  In  fact,  the  system 
performance  for  a  given  inventory  level  can  be  evaluated  exactly  only  for  problems 
having  a  very  small  number  of  different  jobs  and  parts.  System  performance  for  a 
specific  stocking  policy  can  always  be  determined  by  Monte  Carlo  methods,  how¬ 
ever.  When  job  arrivals  form  a  renewal  process  upper  and  lower  bounds  on  system 
performance  can  be  obtained,  which  are  expressed  as  simple  functions  of  the  item 
stock  levels.  The  simple  forms  of  the  bounds  allow  optimization,  subject  to  bud¬ 
get  constraints  and  other  linear  constraints,  such  as  inventory  space.  Optimization 
using  bounds  is  important  for  this  inventory  system  because  it  is  observed  that 
objective  functions  involving  the  exact  performance  criteria  are  neither  convex  or 
concave  in  the  stock  levels.  Several  example  calculations  are  included  in  the  paper 
to  illustrate  the  relationships  between  the  bounds  and  more  precise  calculations. 


2  The  Model 


2.1  Specifications 

The  key  demand  data  for  problems  of  the  type  we  are  considering  are  contained 
in  the  "job  matrix”  J,  which  is  defined  as  follows: 

Jij  =  the  number  of  items  of  type  t  required 

to  satisfy  a  request  (complete  a  job)  of  type  j. 

[In  production  applications,  the  concept  of  a  repair  job  is  replaced  by  a  final 
assembled  product  that  contains  some  collection  of  items  or  subassemblies.] 
Assumption  1: 

Jobs  or  requests  for  items  are  assumed  to  have  random  independent  interarrival 
times  with  mean  1/A.  The  marginal  probabilities  for  job  types  are  defined  as 
follows: 

Pj  =  P{an  arriving  job  or  request  is  of  type  /}. 

The  type  of  each  job  is  assumed  to  be  independent  of  the  types  of  past  and  future 
jobs.  The  marginal  demand  rates  for  individual  items  A,-  t  =  are  then 

determined  by: 


a,=a5:j,,p,.  (1) 

i 

where  Aj  =  average  demand  rate  per  unit  time  for  item  t,  t  =  1, . . . ,n. 

For  any  fixed  time  t,  we  define  the  random  variables 

=  number  of  requests  (jobs)  of  type  j  arriving  by  time  t,  j  =  1, . . .  ,m 

=  number  of  items  of  type  t  demanded  by  time  t,  t  =  1, . . .  ,n. 

Letting 

X{t)  =(Xi(t),...,Ar„(0) 

N{t) 

we  clearly  have  the  matrix  relationship 

X(l)  =  (2) 
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Let  Af(t),  t  >  0  denote  the  counting  process  of  job  arrivals.  That  is,  M(t)  Is 
the  number  of  jobs  which  have  arrived  by  time  t;  i.e., 

Throughout  this  paper  we  assume  M(t)  is  a  renewad  process,  which  restarts  at  each 
restocking  point.  The  J  matrix  introduces  dependence  among  the  components  of 
X'(t).  Even  when  each  job  uses  a  single  distinct  part,  the  components  of  X(t)  may 
be  correlated.  We  summarize  the  dependence  among  the  components  of  N(t)  and 
X(t): 

For  any  t  >  0 


Cov(N,(t),Nj(t))  =  p,pj(Var(M(t))  - 
If  Var(Af(t))  —  E(M(t))  >  0  then  for  all  i,j 

Cov(Xi(t),Xj(t))  >  0. 

2.2  Defining  Performance  Measures 

In  this  subsection,  we  derive  expressions  for  performance  measures  as  a  function 
of  the  initial  stock  level.  These  are  defined  in  terms  of  two  random  variables: 

r  =  time  of  the  first  stockout 
<7  =  number  of  jobs  completed  (demands  met) 
before  stockout 

The  performance  measures  are: 

•  probability  first  stockout  is  after  time  t  P(r  >  t) 

•  expected  time  to  first  stockout  £(r) 

•  probability  that  at  least  k  jobs  completed  before  stockout  P{<7  >  k) 

•  expected  number  of  jobs  completed  before  stockorit  E((t) 

The  most  appropriate  measure  depends  upon  the  application.  For  mission  comple¬ 
tion  problems,  P(r  >  t}  would  be  appropriate  for  a  mission  of  given  duration  t,  or 
E{t)  if  the  duration  is  uncert^.  For  field  repur  applications  E{a)  or  P{a  >  k} 

S 

k 


would  be  appropriate,  because  these  determine  the  probability  that  a  repairman 
can  complete  a  given  number  of  calls  without  stocking  out.  In  the  field  repair  set¬ 
ting,  the  occurrence  of  a  stockout  for  a  single  part  forces  the  repairman  to  return 
to  the  warehouse  for  the  needed  part  and  to  restock  his  kit  at  the  same  time.  For 
the  problem  of  repairing  ships  at  sea,  a  stockout  in  a  single  critical  part  entails 
a  return  to  port  or  some  form  of  emergency  delivery.  For  specified  item  holding 
costs  hi,  these  performance  criteria  lead  to  corresponding  stock  level  optimization 
problems.  Examples  include  mininuzing  the  inventory  cost  of  achieving  a  specified 
performance  level,  or  maximizing  the  performance  that  can  be  achieved  with  an 
inventory  budget  constraint  and/or  space  construnts. 

For  a  given  vector  e  =  ei,...,e„  of  stock  levels,  where  s,-  =  initial  stock  level 
for  item  type  t,  r  the  time  to  first  stockout  is  defined  by: 

r  =  inf {t :  X,(t)  >  «<  +  1  for  some  i,  i  =  1, . . . ,  n}. 

In  words,  r  is  the  first  time  that  the  demand  for  some  part  or  parts  strictly  exceeds 
supply.  Thus 


P{A11  requests  filled  up  to  time  0  =  P{Ar(t)  <  s}  =  P{r  >  t}.  (3) 

When  the  time  horizon  is  not  known  with  certainty,  a  reasonable  performance 
criterion  is  the  expected  time  to  first  stockout,  E(t), 

E{t)=  rP{T>t}di. 

Jq 

Wald’s  equation  provides  a  relationship  between  E{t)  and  E{c).  Define  a  by 

c  =  min{il; :  Xi{k)  >  s,-  +  1  some  1  <  *  <  n} 

Let  Y  denote  the  (random)  time  between  arrivals  of  jobs,  where  E{y)  =  1/A. 
Wald’s  equation  yields  the  identity 

£{r)  =  E{Y)E(c)  =  E{o)/X.  (4) 


Ut 


K,W  = 


number  of  jobs  (demands)  of  type  j  in  the 
first  k  arriving  jobs  since  restocking. 


The  marginal  probability  distribution  for  Nj{k)  is  binomial 

i>{W,(l)  =  «}=  I 

Let  Tk  denote  the  time  of  arrival  of  the  job.  We  denote  the  functions  of  the 
discrete  index  k  by  Nf{k)  =  iV,(r*)  and  X<(ife)  =  X,(rk). 

For  fixed  k,  the  distribution  of  jobs  by  type  is  a  multinomial,  since  the  type  of 
an  arriving  job  is  independent  of  both  its  interarrival  time  and  its  neighboring  job 
types.  Thus, 

/-{,>*}=  Y.  (5) 

»€/*nJr  *»••••  *"»• 

where  /*  =  {s|  E  *i  =  *,  a-  >  0}  and  H  =  {*|s  >0tJz<  s}.  Now, 

i 


k=0M€Hnli 


^1"  •  •  •  Zfn‘ 


Equation  (5)  can  be  used  to  obtain  P{t  >  t}  via  the  identity. 


where  the  last  equality  follows  since  the  job  types  are  independent  of  the  job  arrival 
process. 

Equations  (5)  and  (6)  can,  in  principle,  be  used  to  calculate  the  performance 
criteria  for  a  given  stock  level.  The  calculations,  however,  are  quite  time  consuming 
even  for  fairly  small  numbers  of  items  and  stock  levels.  The  difficulty  is  compounded 
when  we  seek  optimal  stock  levels,  a  process  which  nu^  require  evaluating  a  large 
number  of  policies.  Furthermore,  computational  experience  has  demonstrated  that 
E{a)  and  P{c  >  1}  are  not  concave  functions  of  the  stock  level  s.  Thus  global 
optimality  cannot  be  guaranteed,  even  when  local  optima  are  located. 


2.3  Bounds 

In  this  section,  we  develop  a  series  of  upper  and  lower  bounds  for  the  per¬ 
formance  measures.  These  bounds  are  important  for  two  reasons:  (1)  the  exact 
peffcmnance  measures  are  tedious  to  evaluate  and  for  large  problems  the  computa¬ 
tion  becomes  infeasible  (2)  the  exact  performance  measures  do  not  produce  concave 
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or  convex  optimization  problems,  respectively,  when  maximizing  performance  sub¬ 
ject  to  a  budget  constraint  or  minimizing  cost  subject  to  a  performance  constraint. 
The  bounds  we  devebp,  however,  do  lead  to  concave  or  convex  problems,  as  a  func¬ 
tion  of  the  stock  levels.  Thus  the  bounds  can  play  the  role  of  surrogate  objective 
functions  that  allow  approximately  optimal  stock  levels  to  be  easily  determined. 
Furthermore,  when  the  bounds  are  known  to  lie  below  the  true  performance  level, 
the  optimized  bound  can  be  set  so  that  a  performance  target  is  guaranteed  to  be 
met.  When  several  such  lower  bounds  can  be  computed,  as  shown  in  later  examples, 
the  best  (greatest)  lower  bound  can  be  selected  as  the  best  estimate  of  performance. 

The  power  of  the  bounds  is  increased  in  stock  level  optimization  problems  by 
combining  them  with  more  accurate  (and  time  consuming)  Monte  Carlo  evaluations 
of  the  exact  performance  measures.  That  is,  we  use  the  bounds  in  defining  well 
behaved  objective  functions  for  determining  “good”  stock  levels.  We  then  use  Monte 
Carol  to  evaluate  the  exact  criterion  P{o  >  t},  E{a)  etc.  at  that  stock  level. 
Based  on  this  evaluation,  the  constraints  in  the  surrogate  optimization  could  then 
be  adjusted  and  a  new  improved  stock  determined  that  more  closely  matches  the 
exact  performance  level  target. 

We  begin  with  upper  and  lower  bounds  on  expected  time  to  stockout. 

Lemma  1  With  only  A$»umption  1,  we  have  the  following  upper  bound  on  expected 
time  to  ttoekout 


where 


Proof;  Define 


<T<  =  inf{lfe  :  Xii(ib)  >  s<  -h  1}. 


Clearly, 


Eo  =  E  minfoi, . . . , 

It  follows  from  Jensen’s  inequality  that 

Eo  <  inin{£oi, . . . ,  Ecn). 


(7) 
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The  process  Xi(^)  experiences  an  increment  of  Jii  with  probability  pj.  If  we  let 
Zi,.. .  denote  the  sequence  of  increments  then  we  have 


«i  +  1  <  <  a<  +  1  +  max{7„}, 


where  the  term  on  the  far  right  is  the  maximum  amount  by  which  X,(<7)  will  exceed 
+ 1.  Taking  expectations  on  both  sides  and  applying  Wald’s  equation  again  yields 

mi  -  -  — w:)  • 

m 

and  E[Zi)  =  J,,.  This  yields 

i=i 


Q.E.D. 


The  bound  in  (8)  is  simple  to  compute.  Thus  the  right  hand  side  of  (7)  is  easily 
determined. 

It  is  often  more  important  to  obtain  a  lower  bound  on  time  to  stockout.  The 
lower  bound  provides  a  ’’pessimistic”  estimate  of  performance.  This  is  useful  be¬ 
cause,  by  purchasing  sufficient  inventory,  the  pessimistic  bound  can  usually  be 
brought  to  any  specified  level  of  system  performance.  This  assures  that  acceptable 
performance  will  be  achieved,  although  some  excess  inventory  investment  is  likely 
because  of  the  approximate  nature  of  the  bound. 


Lemma  9  With  prohaiUity  1,  a  tht  number  of  joht  until  etoekout  eatiefiea 

c,  <  o. 


inhere 

and  Rj  denote*  the  maximum  number  of  eorueeutive  job*  of  type  j  which  can  be 
*erved  before  atoekout, 

Ej  =  mm  ; — ). 


0 


Prc  i:  It  is  sufficient  to  show, 


{w  :  JN{k)  <  «  +  1}  2  <  !}• 

i=i  "i 


Suppose  that  for  some  u, 


<  1,  the  definition  of  Rj  then  gives 
i=i 


Now 


Thus 


or 


V _ - _ 


-NAk)  <  1 


(  min  =  max{  } 


js:l  «•  T 

fl 

52  Ji}Nj{k)  <  s,  +  1.  t  =  1, . . . ,  n, 
»=i 


Q.E.D. 


Lemma  2  produces  a  major  simplification  because  the  vector  relationship  JN{k)  < 
«  +  1  is  reduced  to  a  one  dimensional  relationship  expressed  as  a  weighted  sum  of 
binomial  random  variables.  Lemma  2  obtains  a  pessimistic  bound  for  the  proba¬ 
bility  distribution  of  the  number  jobs  that  are  completed  before  stockout.  In  some 
applications,  the  probability  distribution  of  the  time  until  stockout  is  a  more  appro¬ 
priate  performance  measure.  It  can  be  noted  in  the  proof  of  Lemma  2,  that  all  the 
implications  are  equally  valid  if  Nj{k)  is  replaced  with  the  number  of  jobs 

of  type  j  arriving  by  time  t.  Thus,  using  arguments  analogous  to  those  in  Lemma 
2,  we  can  establish  the  following  lower  bound  on  the  probability  distribution  of  the 
time  to  stockout. 


Lemnui  S 


where 


and 


<•}  =  P{r  >  t}  >  P{r,  >  t} 
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The  situation  in  Lemmas  2  and  3  is  depicted  graphically  in  Figure  1.  Suppose 
there  are  2  possible  jobs  and  2  parts.  The  first  type  of  job  involves  part  1  only  and 
the  second  parts  1  and  2.  Suppose  we  stock  3  of  part  1  and  2  of  part  2.  Plotted 
along  each  axis  are  the  number  of  jobs  of  type  1  and  2  observed.  The  solid  lines 
depict  the  constraints  Jx  <  s  +  1.  Thus  the  lattice  points  within  this  region  indicate 
combinations  of  jobs  which  can  be  repaired  from  stock  on  hand.  The  dashed  line 

m 

is  the  set  of  points  such  that  ~  process  N{k)  moves  from  lattice 

*•=1 

point  to  lattice  point  in  a  random  fashion,  but  always  to  the  northeast.  (A  typical 
realization  is  given  by  the  line  marked  with  arrows).  The  lemma  states  that  the 
N{t)  process  always  crosses  the  dashed  lines  before  crossing  the  solid  boundary. 
Lemma  2  thus  implies  that  the  complementary  cumulative  distribution  of  o.  lies 
to  the  left  of  that  of  a  in  the  sense  that  P{a,  >  k}  <  P{a  >  k}  all  k.  Thus  the 
distribution  of  a.  provides  a  bound  on  the  distribution  of  a. 


2.4  Computing  the  Bounds 


Probability  of  Stockout 

The  distribution  of  a,  is  tedious  to  calculate  exactly,  but  it  is  considerably  easier 
to  compute  than  the  distribution  of  c.  Define 


Then 


Pr{0,  >k}  = 


i=i 

E-  -E 


kl 


-pST- 


(9) 


KHi  t\\  , , .  fnt' 

For  large  values  of  k,  {P{a,  >  k}  can  be  approximated  by  the  tail  of  a  Normal 
distribution.  Since  the  Nj(k)  j  =  are  multinomially  distributed,  we  let 

,=i  >=i  "j 

»■!  ye«+l 

Approximating  the  multinomiad  by  a  multivariate  normal  with  the  same  mean  and 
covariance  matrix  yields 

P{c.  >  *}  =  ^  S  J) » 
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where  $(•  |a,  6)  represents  a  univariate  normal  distribution  with  mean  a  and  variance 
6. 

The  Normal  approximation  can  also  be  used  to  compute  P{t,  >  t}  = 

<  1}  for  Lemma  3.  Here  we  simply  replace  Nj(k)  above  with 
Nj{t)  and  determine  the  revised  mean  and  variance  n{t)  and  We  then  have 

by  the  Central  Limit  Theorem  that 

This  approximation  is  best  when  t  is  such  that  the  Nj(t)  are  fairly  large.  This  in 
turn  implies  that  the  s,-  should  be  large  enough  to  serve  a  substantial  number  of 
jobs.  Thus,  the  Normal  approximations  are  appropriate  when  both  the  stock  levels 
and  the  number  of  jobs  to  be  completed  before  stockout  are  large. 

Expected  Time  to  Stockout 

The  expected  time  to  stockout  E{t)  can  be  bounded  below  by  using  the  prob¬ 
ability  distribution  P{r.  >  t}.  That  is, 

E(r,)  =  jT*  P{r.  >  t}dt  <  P{t  >  i)dt  =  E(r).  (10) 

Thus  E(t,)  provides  a  lower  bound  on  E{t).  Given  the  complexity  of  the  analytical 
expression  for  the  distribution  of  o,  (and  hence  r.),  the  integral  in  (10)  may  be 
difficult  to  calculate.  We  can,  however,  develop  a  simple,  precise  approximation  to 
E{p,)  (and  hence  to  £(r.)). 

Lemma  4  Tht  pttnmiatie  bound  E{a,)  on  number  of  jobe  until  stockout  eatiefiea 

^  <  (1  +  mw{l/P,>]/(J3p,/Ji,). 

i=i  ' 

Proof;  Define  W*  by 

Wt  =  N{k)-pk 

where  p  =  {pi, . . .  ,Pm);  Wt  is  a  discrete  time,  P"  valued  martingale.  To  see  this, 
note 

E{Wm  I  W*)  =  E{W,  +  WM  -  p  I  W,), 

where  P{u>k+i  =  e<)  =  Pk*  representing  the  s'*  unit  vector  in  It  follows  then 
that 

E{Wm\W^)  =r(W*|W0  +  r(ti;»+i-p) 
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For  each  n,  ||  Wi,  ||<  n  establishing  the  martingale  property.  Note  that  since  N{k) 
moves  only  to  the  northeast  (increases  according  to  the  usual  vector  partial  order) 

fb 

this  implies  that  P{a,  <  oo)  =  T;  hence,  a,  <  ^(si  +  1).  The  optional  sampling 
theorem  applied  to  Wu  yields 


0  =  E{W^.)  =  E{N{a,))  -  pE{a,).  (11) 


Thus  £(N(o,))  =  pE{c,).  It  follows  from  the  definition  of  a„  however,  that  with 
probability  1 

i=i  "i  * 

Taking  inner  products  on  both  sides  of  the  vector  equation  (12)  with  respect  to  the 
vector  (l/i2i, . . . ,  1/Am)  and  applying  (13)  yields 


1  <  ^  1  +  max(l/A,).  (13) 

j=i 

Where  the  last  term  of  (15)  bounds  the  amount  by  which  N{k)  may  overshoot  the 

ft 

boundary  of  the  set  <  1.  Solving  (14)  yields  a  lower  bound  L 


(14) 

Q.E.D. 


Note  that  in  the  extreme  case  in  which  each  job  uses  every  part,  the  bound  L  is 
tight.  In  view  of  equation  (7),  L/X  is  a  lower  bound  for  E(r).  Summarizing,  we  have 
obtained  the  following  bounds  and  apprcodmations  for  the  performance  measures. 
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Table  2.1.  Summary  of  Bounds  Obtained 


£(») 

expected  no  of  jobs  completed 

Lower 

(1  -1-  max,{l/J^})/  Pi/Ri 

Upper 

P{X{t)  <  5}( 

=  P{no  stockout  by  t}|  =  P{t  >  t) 

Lower 

p{zT-im)/Ri  >  1) 

P{a  >  k}|  =  P{at  least  k  jobs  served}] 

Lower 

<  1} 

=  ♦(!/)•(*);,»’(*))  for  Urge 

where: 

—  niin{,7jy>o){(s»  + 

li 

«^(*)  =  *E(i/«,)*P,(l-p,)-2*2:  E  1/Rii/RiPiPi 

t«i  iw+i 

S  Finding  Approximately  Optimal  Stocking  Policies 

l|he  goal  of  our  inventory  analyiis  is  to  provide  "good”  policies  for  stocking.  In 
order  to  do  this,  we  need  to  evaluate  the  performance  of  a  large  number  of  policies. 
This  is  a  difficult  task  in  view  of  the  complexity  of  equations  (5)  and  (6). 

This  section  uses  the  heuristic  of  replacing  the  actual  performance  criterion  by 
our  bounds  on  the  performance  criterion.  The  bounds  therein  provide  surrogate  op¬ 
timisation  problems  whose  solutions  lead  to  Vxxl”  inventory  policies.  The  bounds 
specified  above  offer  particularly  easy  optimisation.  In  addition,  optimisations  us¬ 
ing  the  lower  bounds  offer  guarantees  of  minimum  performance.  In  our  example 
^plkatidn,  we  concentrate  on  the  criterion  of  the  expected  time  until  stockout. 

Consider  the  problem  of  stocking  a  mobile  repairman  who  operates  in  the  fol- 
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lowing  way.  Starting  with  an  initial  kit  of  parts  (ai, . . . ,  an)i  he  answers  repair  calls 
one  at  a  time  until  encountering  a  job  which  cannot  be  repaired  from  the  remaining 
contents  of  his  kit.  At  this  point  he  returns  to  the  parts  depot  to  obtain  the  parts 
for  that  job  as  well  as  to  restock  the  kit.  The  trip  to  the  parts  depot  results  in 
lost  work  time  and  customer  good  will.  We  take  as  our  objective  the  maximisation 
of  the  expected  number  of  repair  jobs  that  can  be  completed  before  restocking  is 
required.  Space  is  often  tightly  construed  in  such  service  systems  since  the  repur- 
man’s  part  supply,  as  well  as  tools,  must  fit  inside  his  car.  In  addition,  budgetary 
constraints  are  an  important  consideration.  Our  objective  will  be  to  maximize  the 
expected  number  of  repair  jobs  completed  between  restockings,  subject  to  the  bud¬ 
get  and  space  constraints.  Let  the  cost  of  part  *  for  inventory  purposes  be  hi  and 
the  space  required  by  n.  Then  the  problem  to  be  solved  can  be  stated  as 

max  £(r)  (15) 

s.t.  ^ 

Em.  <v 

’n* 

En*.  $  -R- 

fisl 

This  example,  as  well  as  similar  practical  problenos,  can  be  written  as 

max  ^  E{t)  (16) 

S.t. 

As  <b 
s  >0, 

where  A  is  a  matrix  containing  the  coefficients  of  the  constraints  on  the  stock  levels. 

In  (16),  E{t)  is  calculated  from  E{t)  =  E{o)fX.  The  previously  discussed  diffi¬ 
culty  in  calculating  E{o)  is  compounded  in  the  problem  of  finding  an  optimal  stock 
level  vector;  any  approach  to  optimization  by  enumeration  requires  the  evaluation 
of  E{a)  for  a  large  number  of  stock  level  vectors.  The  possible  number  of  such 

vectors  grows  roughly  in  proportion  to  JJsj,  where  (si«. •..#«)  are  the  maximal 

feasible  stock  levels.  It  b  interestii^t  to  note  that  the  problem  of  maarimlKltn  the 
expected  time  to  stockout  b  difficult  even  when  the  jobe  consist  of  ffiatinct  groups 
of  parts  (i.e.,  when  part  demands  are  independent).  Thb  b  so  because  the  time  to 
stockout  b  the  minimum  time  to  stockout  ctf  each  of  the  parts. 


15 


Because  of  the  difficulty  encountered  in  repeatedly  evaluating  £(o),  we  propose 
two  heuristic  approaches  to  optimisation.  Our  approaches  involve  optimizing  the 
upper  and  lower  bounds  for  £(o).  Let  V (jj,  denote  the  value  of  the  problem 

(16)  at  an  optimal  policy.  We  define 


=  ^max  U(si . s„)(l/A) 


(17) 


s.t. 

As  <b 
s  >0 

where  U  is  given  in  (8)  (the  dependence  on  Si,...,Sn  is  made  explicit  in  (17)  for 
clarity).  Rewriting  (17)  yields 


max  {min{ 
(*i . *.)  i 


Sj  -f  max,,{J,j)  +  1 


mm 


»=i 


As  >  6, 

Si  >  0. 

This  expression  for  can  be  rewritten  by  introducing  additional  inequalities: 


s.t. 


max  t/A 


(18) 


Y,  ^  max<,{  J<,}  +  1 

As  <b 

Si  >  0,t  =  l,...,n. 


R,  as  seems  prudent  pven  the  approximate  nature  of  (18),  we  ignore  the  integer 
restrictions  on  the  s<,  (18)  is  a  standard  linear  program.  A  moment’s  reflection 
reveals  that  this  is  exactly  the  linear  program  one  would  obtain  by  ignoring  the 
random  nature  of  the  breakdowns  and  aohring  a  deterministic  stocking  problem  in 

m 

which  each  breakdown  required  Y^*iPi  ** 

i-i 

On  the  other  hand,  maximising  our  lower  bound  (equation  (14)),  g^vcs  a  con¬ 
servative  stocking  policy  and  yields  the  optimisation 
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(19) 


=  . . 

s.t.  As  <  b 


Rj  ^  (««  +  l)Mi  c&ch  i  such  that  >  0,  all  j 

*i  >  0i*  = 


which  is  equivalent  to 


(20) 


s.t.  Rj  <  (si  +  1)/Jij  each  i  such  that  Jij  >  0  all  j. 
i4«  <6 

Si  >  0,»  =  l,...,n. 

Problem  (20)  can  be  solved  as  with  linear  constraints  and  a  separable  convex  ob* 
jective. 

Intuitively,  is  the  result  of  optimizing  the  “certunty  equivalent”  problem  in 
which  the  random  part  demands  are  replaced  by  deterministic  demands  equal  to 
the  expected  demand  per  job.  On  the  other  hand  represents  a  weighted  max> 
min  criterion:  the  number  of  each  job  that  can  be  served  is  determined  by  the  part 
that  runs  out  of  stock  soonest.  The  time  to  stockout  is  estimated  by  the  weighted 
average  of  these  stockout  times.  It  seems  reasonable  that  the  policy  given  by  the 
optimization  of  the  lower  bound  will  be  good  when  there  is  a  high  level  of  part 
commonality  between  jobs.  Indeed,  in  the  extreme  case,  in  which  each  part  is  iised 
in  every  job,  gives  the  optimal  stocking  policy.  For  the  other  performance  criteria 
in  Table  2.1,  conservative  stocking  policies  can  be  obtained  in  a  similar  manner  by 
using  the  lower  boimds  to  form  surrogate  objective  functions. 

4  Examples 

We  now  tom  to  specific  illustrative  numerical  examples  of  the  applications  de- 
scribed  in  section  4.  The  nature  of  the  solution,  as  well  as  the  efficacy  of  our 
approaches  to  it,  depend  critically  on  the  nature  of  the  job  matrix  J.  Our  specific 
examples  are  designed  to  illustrate  typical  situations  which  may  occur  in  practice. 
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rather  than  the  "average”  performance  of  our  techniques  applied  to  "randomly  gen¬ 
erated”  problems.  Thb  approach  will  give  some  qualitative  insights  into  the  effects 
of  various  types  of  part  demand  dependencies. 

Our  numerical  examples  concentrate  on  the  problem  of  maximizing  the  expected 
time  to  stockout,  subject  to  a  budget  constr^t.  We  consider  3  policies,  the  policies 
resulting  from  optimizing  the  upper  and  lower  bounds  and  the  policy  that  results 
from  a  part-fill  type  calculation,  e.g.,  in  which  the  dependence  between  parts  is 
ignored  and  stock  levels  are  set  so  that  the  minimal  expected  time  to  stockout  is 
maximized  for  each  of  the  marginal  demand  processes  (subject  to  a  budget  con¬ 
straint).  This  is  obtained  from  the  optimization 


max  f/A 

s.t.  t<-sr^  *  =  1 . n  (21) 

As<b 

*>0,s,>0  t=;l,...,n 

The  problem  (21)  is  similar  to  the  upper  bound  optimization  (18)  (the  righthand 
sides  are  changed  by  a  constant).  The  solution  to  problem  (21)  has  the  effect  of 
tn>-riTni«ing  the  tniTtitmim  expected  time  to  achieve  a  0  stock  level  for  each  of  the 
marginal  demand  processes  considered  separately.  As  it  concentrates  only  on  the 
marginal  part  demands  we  call  this  the  "part-fill”  heuristic. 

In  order  to  minimize  the  integer  effects  caused  by  part  usages  which  occur  in 
batches,  each  of  the  examples  the  part  usage  matrix  J  was  chosen  to  be  a  0-1 
matrix.  For  simplicity,  we  assume  that  the  matrix  A  consists  of  a  single  budget 
constrdnt  of  ^  hi»i  <  b  where  h,  are  the  unit  holding  costs. 

Each  example  has  20  parts  and  10  jobs,  the  job  probabilities  were  assumed 
equal  and  part  costs  for  examples  are  given  in  table  1.  The  job  matrices  are  given 
In  Tables  2  -  7.  In  problem  1  each  job  omsists  cS  two  distinct  parts.  Problem 
2  entdls  a  low  level  of  commonality  between  the  jobs  and  problem  3  a  high  level 
of  commonality.  Problem  4  involves  a  cascading  part  usage  structure,  each  job 
requiring  a  subset  of  the  parts  used  in  previous  (lower  nuinbered)  jobs.  Problem 
6  consists  of  a  mixture  of  the  overlapping  and  non-overlapping  cases.  Parts  1  - 
6  are  common  to  many  jobs;  7  -  20  are  used  in  diffoent  jobs.  Problem  6  has 
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the  same  part  usage  matrix  as  problem  5,  but  the  part  costs  are  changed  so  that 
hi  =  . . .  =  =  5  and  hy  &=  h|  «  . . .  =  hjo  =  0.01.  In  each  case  we  take  b  =  500 

and  A  =  1.  Clearly,  our  choices  of  structures  for  J  do  not  exhaust  all  possible 
commonality  types.  They  are  meant  to  represent  in  a  general  way  the  structure  of 
part  commonalities  experienced  in  practice.  Such  structures  can  be  characterized 
as  “separate”  (problem  1),  “minimal  overlap”  (problem  2),  “significant  overlap” 
(problem  3),  “cascading”  (problem  4),  and  a  mix  of  “overlapped”  and  “separate” 
(problem  5  and  6}. 

Table  7  compares  the  mean  and  variance  of  the  time  until  stockout  obtained 
via  Monte  Carlo  simulation.  In  each  case  we  took  the  time  between  jobs  to  be 
exactly  1  (A  =  1].  The  policy  obtained  using  the  upper  bound  as  a  surrogate 
objective  fimction  was  generally  dominated  by  the  policy  obtained  from  the  “part- 
fill”  optimization.  The  policy  obtuned  from  optimizing  the  lower  bound  performed 
relatively  well  in  those  cases  in  which  there  was  a  high  level  of  conunonality  of 
parts  between  the  jobs,  and  relatively  poorly  when  there  was  little  commonality. 
While  increasing  the  level  of  part  commonality  between  jobs  improved  the  relative 
performance  of  the  lower  bound  heuristic,  part  commonality  alone  is  not  enough  to 
cause  the  lower  bound  heuristic  to  perform  much  better  than  the  part  fill  heuristic. 

Problem  6  demonstrates  the  conditions  under  which  the  lower  bound  heuristic 
performs  better  than  either  the  upper  bound  or  the  part-fill  heuristic.  In  problem 
6,  those  parts  common  to  several  jobs  (l  -  6)  are  expensive  and  those  parts  used 
in  a  single  job  (7  -  20)  are  cheap.  No  job  can  be  performed  without  at  least  one 
part  specific  to  that  job  alone.  By  concentrating  on  the  marginal  demand  rates,  the 
upper  bound  and  part-fill  heuristics  over  invest  in  those  parts  common  to  several 
jobs  and  under  invest  in  the  unique  parts.  In  contrast,  the  lower  bound  heuristic 
equalizes  the  number  of  each  type  of  job  to  be  done  and  thus  equalizes  its  parts 
expenditure.  For  each  of  the  problems  tested  the  lower  bound  heuristic  resulted  in 
nearly  identical  part  stocking  levels  across  parts,  and  the  upper  bound  and  part- 
fill  heuristics  resulted  in  widely  varying  stock  levels.  This  is  due  in  part  to  the 
assumed  equal  job  arrival  rates  and  in  part  to  the  “max-min”  nature  of  the  lower 
bound  heuristic. 

To  further  test  the  hypothesis  that  the  quality  of  the  solution  obt^ed  from 
the  lower  bound  improves  as  the  problem  becomes  more  dependent,  we  generated 
2  sets  of  10  problems  with  randomly  constructed  part  usage  matrices.  The  first  set 
had  an  average  of  2  parts  per  job,  the  second,  12  parts  per  job.  Table  8  summarizes 
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the  results  of  the  experiment.  Once  again,  the  solution  based  on  using  the  lower 
bound  as  a  surrogate  objective  fared  better  when  there  was  a  high  level  of  part 
commonality  between  jobs.  We  also  see  that  the  variance  of  the  time  to  stockout 
was  smaller  for  the  lower  bound  heuristic  policy  them  for  the  part-fill  policy.  The 
reason  for  this  stems  from  the  structure  of  the  objective  function  in  (21).  With  a 
high  level  of  part  commonality,  problem  (21)  tended  to  produce  policies  for  which 
the  number  of  each  type  of  job  which  could  be  completed  from  the  kit  is  similar. 
For  a  policy  in  which  the  number  of  each  kind  of  job  completable  from  the  kit 
is  exactly  the  same,  the  variance  in  time  to  stockout  is  0.  Finally,  it  should  be 
noted  that,  as  expected,  the  quality  of  the  lower  bound  as  an  estimate  of  system 
performance  improved  dramatically  with  higher  levels  of  part  commonality. 

5  Conclusion 

A  key  purpose  of  this  paper  has  been  to  call  attention  to  the  importance  of  the 
dependence  among  item  demands  in  a  multi-item  inventory  system.  The  type  of 
dependence  considered  in  this  paper  arises  naturally  in  many  multi-item  settings, 
as  discussed  in  the  introduction. 

When  part  or  item  demands  show  significant  dependence,  the  potential  error 
arising  from  ignoring  this  dependence  can  be  large.  While  ignoring  the  dependence 
between  item  demands  may  lead  to  serious  errors  in  assessing  system  performance, 
exact  stock  level  optimisation  with  item  demand  dependence  seems  fraught  with 
difficulty.  Evaluation  of  the  exact  formulas  for  system  performance  for  a  single 
stocking  policy  is  computationally  demanding,  even  for  systems  with  very  modest 
size.  Monte  Carlo  is  effective  for  evaluation  of  individual  stocking  policies,  but  is 
not  well  suited  for  optimization  because  the  objective  functions  are  not  concave  or 
convex. 

Our  solution  to  this  problem  is  to  develop  easily  computable  convex  or  concave 
bounds  on  system  performance.  In  Section  4,  these  bounds  were  used  as  siirrogate 
objective  functions  to  find  policies  which  maximized  the  performance  of  the  bound, 
subject  to  a  budget  and  other  constridnts.  The  simplicity  of  the  bounds  gives  rise 
to  relatively  simple  surrogate  optimization  problems.  Monte  Carlo  is  then  used  to 
test  the  actual  performance  of  the  stock  levels  determined  by  true  (q>timisation. 

The  results  of  our  limited  experience  with  the  bounds  are  encouraging.  The 
part-fill  heuristic  gave  reasonable  performance  when  there  was  a  low  level  of  part 


demand  dependence  and  the  lower  bound  heuristic  gives  very  good  performance  for 
high  levels  of  item  demand  dependence.  In  each  case,  the  bounds  themselves  may 
give  rather  large  over-  or  underestimates  of  system  performance,  but  the  policies 
resulting  from  the  optimization  of  the  bound  performed  well  in  the  Monte  Carlo 
experiment. 

Clearly,  many  interesting  and  unsolved  problems  remain  in  analyzing  stocking 
policies  for  interdependent  multi-item  inventory  systems.  Given  the  increased  in¬ 
terest  in  multi-item  inventories  in  production,  as  well  as  in  repair  kit  inventory 
problems,  we  hope  that  the  results  obtained  in  this  paper  will  stimulate  further 
progress  in  this  research  area. 
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Lower  Bound 
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Table  8. 

Montecarlo  Results  Sample  Mean  and  Variance 
of  Time  to  Stockout  -  Randomly  Generated  Problems. 


Lower  Bound 
Mean 

Heuristic 

Variance 

Fill-Rate 

Mean 

Heuristic 

Variance 
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job  (average). 
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44.78 

12  parts  per 
job  (average). 
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